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Abstract. We examine the possibility to extract information about the DN and DN interactions from 
the pd^D'^ D~p reaction. We utilize the notion that the open-charm mesons are first produced in the 
annihilation of the antiproton on one nucleon in the deuteron and subsequently rescatter on the other (the 
spectator) nucleon. The latter process is then exploited for investigating the DA'' and DN interactions. 
We study different methods for isolating the contributions from the D'^p and D~p rescattering terms. 
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. 1 Introduction 

The distortion of charm in nuclear matter remains an 
' heavily discussed issue since the first proposals [1] to use 
, charmonia and open charm as a probe of the early stage of 

■ heavy-ion collisions, for a review see [2] . It was expected [3] 
I that charmed final-state interactions (FSI) either at the 
. partonic or the hadronic rescattering level would not dis- 
' tort the spectra initially produced in heavy-ion collisions, 

because the cross sections for any such (elastic and inelas- 
' tic) scattering processes are sufficiently small. Further- 

■ more, gluon radiation or bremsstrahlung [4-6], which 
distorts the original charm spectrum as well, becomes the 
dominant energy loss mechanism only if the heavy charmed 
quarks are ultra-relativistic. That is similar to the brems- 
strahlung losses of electrons passing through a hydrogen 
target [7,8]. However, in the present experiments a large 
fraction of the heavy quarks are produced with momenta 
less than their mass and therefore the radiation losses 
might be negligible. In that case the heavy charmed quarks 
and antiquarks and, specifically, the finally detected D 
and D mesons are presumably not distorted in the nu- 
clear environment and thus can probe the initial stage of 
the interaction, possibly, the Quark Gluon Plasma. 

On the other hand it was argued [8-11] that the two 
basic processes involved in the energy-loss mechanism of 
charmed particles moving in nuclear matter, namely gluon 
radiation and elastic scattering, might be equally impor- 
tant and non-negligible. Only recently the situation chan- 
ged due to the PHENIX and STAR experiments [12-15] 
at the Relativistic Heavy Ion Collider (RHIC). These new 
measurements indicate a substantial suppression of the 
production of open-charm mesons with transverse mo- 
menta above ~1 GeV/c from central Au+Au collisions, 



as compared to that from d+Au collisions. This observa- 
tion could not be assigned to gluon radiation of charm and 
thus indirectly points to an importance of distortions due 
to the FSI. 

Apparently, the interactions of charmed quarks or an- 
tiquarks in nuclear matter are not the same as the inter- 
actions of open-charm mesons. However, it is clear that 
a reasonable understanding of elastic scattering involving 
particles with charm on a hadronic level is highly impor- 
tant. While one could not measure directly the interaction 
of the charmed and light quarks and antiquarks, the DN 
and DN interactions can be studied experimentally. These 
could serve as a basis to construct phenomenologically the 
charmed FSI at the partonic rescattering level. 

The basic problem is the complete lack of relevant ex- 
perimental data. This situation is expected to change with 
the operation of the Facility for Antiproton and Ion Re- 
search (FAIR) at Darmstadt (Germany). The Proton AN- 
tiproton at DArmstadt (PANDA) Collaboration [16] in- 
tends to investigate the distortion of open-charm mesons 
[17-19] in matter and in the vacuum. The matter mea- 
surements are based on D and D meson production in an- 
tiproton annihilation on different nuclei in order to study 
the A-dependence. The vacuum measurements explore the 
production of open-charm mesons by annihilating antipro- 
tons on the deuteron and, through the rescattering of the 
produced D and D mesons on the spectator nucleon, the 
interaction in the DN and DN systems. 

In the present paper we examine the possibility to 
extract information about the DN and DN interactions 
from the pd-^D^D'p reaction. The study is based on the 
notion that those open-charm mesons are first produced 
by annihilating the antiproton on one of the nucleons in 
the deuteron and subsequently rescatter on the other (the 
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spectator) nucleon. The latter_ process is then exploited 
for investigating the DN and DN interactions. 

To explore the potential of pertinent experiments we 
perform concrete calculations taking into account the nu- 
cleon exchange Born diagram, corresponding to the ele- 
mentary NN — > DD annihilation process, cf. Fig. la), 
as well as rescattering diagrams involving the DN and 
the DN interactions, see Fig. lb). For those interactions 
wc employ realistic DN^DN [20] and DN-^DN [21] 
scattering amplitudes. This is a substantial improvement 
as compared to previous studies [22-24] which relied on 
rather simple assumptions as far as the DN and DN scat- 
tering amplitudes were concerned. 

The paper is structured in the following way: In the 
subsequent section we introduce briefly the formalism used 
for calculating the reaction pd^D^D~p. The utilized in- 
teraction models in the DN and DN channels arc intro- 
duced and discussed in Sect. 3. In particular, we present 
results for total and differential cross sections for the var- 
ious charge channels. In Sect. 4 a short overview on our 
present knowledge on the elementary NN DD reac- 
tion is given. Our results for the reaction pd^D^D~p 
are shown in Sect. 5. Here the emphasis is put on the 
exploration of different methods for detecting and isolat- 
ing the contributions from the D^p and D^p rescattering 
terms. For that purpose we consider the spectator mo- 
mentum distribution, Dalitz plots, the missing mass of 
the exchanged meson and correlations between properly 
defined scattering planes. The paper ends with a brief 
summary. As a test of our approach we also apply it to 
multipion production in pd annihilation and we compare 
our results with data available for the pd— >7r+27r~p and 
pd^2n'^3'K~p reactions. Those results are included in an 
appendix. 

2 Formalism 

In this section we introduce the formalism we use to in- 
vestigate the effects of the DA'' and DN interactions in 
antiproton annihilation on the deuteron. 

Fig. 1 illustrates processes contributing to the reac- 
tion pd^DDN . The diagrams of interest arc the nucleon 
exchange Born diagram, Fig. la), and the meson rescat- 
tering diagram, Fig. lb). The corresponding amplitudes 
will be denoted below by Ta and Tb, respectively. 

In what follows we assume for convenience that the 
spectator nucleon - the nucleon that does not enter the 
annihilation vertices - is the proton. The case of a neu- 
tron as spectator can be treated in a similar manner. The 
nucleon exchange Born diagram. Fig. la), leads to the well- 
known reaction amplitude [25-27] 

Ta ^ tpd{Ps)TA, (1) 

where ipd{Ps) is the momentum-space deuteron wave func- 
tion, with Ps being the proton spectator momentum, and 
Ta is the amplitude for the pn D~D'^ annihilation pro- 
cess. After summation over spin states the squared ampli- 
tude is given as 

\Ta\^ = [u{psf+w{psf]\TAf, (2) 




Fig. 1. Contributions to the reaction pd^DDN: a) The Born 
(nucelon exchange) diagram. Ta denotes the annihilation am- 
plitude, b) Meson rescattering diagram. Tm denotes the niesori- 
nucleon scattering amplitude. Note that both DN and DN 
scatterings contribute to the reaction amplitude. 

where u and w stand for the s- and c?-wave components 

of the deuteron wave function. The overall size of the 
nucleon exchange contribution in Fig. la) is determined 
by the annihilation amplitude Ta, which depends on the 
meson channels produced in the pn annihilation. On the 
other hand, the spectator momentum distribution is gov- 
erned predominantly by the deuteron wave function, as 
indicated in Eq. (2), provided that Ta is a slowly varying 
function of the energy. 

Next we consider the rescattering diagram of Fig. lb) 
where one of the mesons produced at the annihilation 
vertex is scattered off the spectator nucleon. In general, 
the intermediate meson in this diagram is not necessar- 
ily the same as the final rescattered meson. It could be 
an intermediate D* vector-meson, for example. Thus, in 
principle a sum over all possible intermediate states is re- 
quired. But in the following let us regard explicitly the 
features of the rescattering mechanism involving elastic 
Mp^Mp scattering only. Note that the formalism can 
be easily extended to other, non-diagonal transitions. We 
average over the spins in the annihilation and scattering 
vertices and take into account only the s-wave component 
u{ps) of the deuteron wave fimction. The d-wave compo- 
nent is expected to play a much less important role for the 
rescattering contribution [27, 28] and, therefore, we ignore 
it here in this exploratory study. The integration for the 
rescattering diagram runs over the three momentum of the 
spectator proton in the loop and can be split into on-shell 
and off-shell parts, which we denote by T^" and T°^^. The 
on-shell part is defined by taking the intermediate proton 
to be on-shell and is given as [27,29-32] 

9+ 




dqu{q)^TMTA, (3) 
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where Tm is the meson-proton scattering anipUtude, q is 
the internal proton loop momentum, Eg = (q^ + m^)^/^ 
with rrip the proton mass, and p is the sum of the momenta 
of the final proton and the rescattered meson. The limits 
of the integral are given by [33] 



Q± = 



\P\ 



E* ± 



E 

"T72 



(4) 



where E=Ei\i+Ep is the sum of the energies of the final 
proton and the rescattered meson, while s]^^ = {E'^—p^y/'^ 
is their invariant energy and 
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(5) 



with thm the meson mass. The evaluation of Ti, requires 
the knowledge of both amplitudes Tm and Ta within the 
range of to allowed for the considered reaction. 

The effect of the off-shell part of the rescattering inte- 
gral was investigated in detail in Refs. [27,29-31,34,35]. 
In those studies it was found that the shape of the specta- 
tor momentum distribution, as given by the on-shell part, 
remains essentially unchanged when the off-shell contri- 
bution is added. At the same time, the magnitude of the 
off-shell contribution is significant and can lead to modifi- 
cations of the on-shell results in the order of 30% or more, 
but depends strongly on the specific off-shell behaviour of 
the annihilation and scattering amplitudes. Explicitly, the 
off-shell part of the amplitude can be written as [27,31, 
34] 



rpoff 



00 

— ^ / dq u{q) Tm Ta 
327r^p J Eq 



X In 



EM+Ep—Eq—E^ 



EM+Ep — Eq — Ej^ 



(6) 



with 



any really significant enhancement in that spectator dis- 
tribution with respect to the results based on the Born 
term can be definitely attributed to meson rescattering 
processes on the spectator. Since the scattering amplitude 
T}) is directly proportional to the Mp—^Mp amplitude Tm 
and the annihilation amplitude Ta enters both the Born 
and rescattering diagrams, the enhancement in the high 
momentum part of the spectator distribution is related to 
the relative magnitudes of their contributions and, thus, 
essentially driven by the scattering amplitude Tm- 

In any case, the isolation of the effects of Tm from 
the data is by no means a trivial task. There are also 
uncertainties related to possible contributions from other 
processes besides those considered above, as well as uncer- 
tainties related to the evaluation of off-shell corrections to 
the scattering amplitude beyond Eq. (6). Nonetheless, the 
situation is still manageable since those uncertainties can 
be kept under control to some extent by selecting carefully 
the reaction kinematics - see for example Ref. [36] . 

Our predictions for the pd^DDN reaction are pre- 
sented and thoroughly discussed in Section 5. In the Ap- 
pendix we illustrate the applicability of the formalism 
to the proton spectator distributions as measured in the 
pd^'!T^2'K~p and pd-^27:^ 3tt~ p reactions. Besides provid- 
ing support for the reliability of the approach adopted for 
the present investigation, those results are also useful for 
revealing similarities but also differences in the utilization 
of the reaction pd nM N for exploring the MN inter- 
action for different mesons and/or different kinematical 
regimes. 

All calculations are done with the deuteron wave func- 
tion of the CD-Bonn NN potential [37]. But for shedding 
light on the model dependence we also emply the wave 
functions of the Paris [38] and full Bonn [39] NN po- 
tentials. We ignore the effects of the Coulomb interaction 
in this exploratory study. Furthermore, we use averaged 
masses. Specifically ,_ we use niD = 1866.9 MeV for the 
mass of the D and D mesons. 



3 The DN and DN interactions 



2 , 2 



(7) 



where m,px is the mass of the exchanged meson. In the 
present exploratory study we do not take into account 
those off-shell contributions. We stress again that these 
are highly model-dependent in any case. 

Let us emphasize the general features of the reaction. 
It is clear that the spectator proton momentum spectrum 
is sensitive to the reaction mechanism. The nucleon ex- 
change Born term |Ta|^, given by Eq. (2), produces low 
and high momentum components, where the high momen- 
tum part of the spectator distribution is dominated by 
w(gr), the d-wave component of the deuteron wave func- 
tion. The d-wave component is more strongly model de- 
pendent and, therefore, the corresponding contribution of 
the Born term to the high momentum part of the spectator 
distribution is affiicted with some uncertainties. However, 



Considering charm and charge conservation antiproton- 
deuteron annihilation allows to study the following D- 
meson production reactions 



pd D°D°n, 
pd D+D-n, 
pd D°D-p . 



(8) 

(9) 
(10) 



They involve the DN as well as the DN scattering ampli- 
tude. Here we investigate only the pd-^D^D^p reaction. 
We select this channel because, in principle, all final parti- 
cles can be measured for this reaction, which implies that 
there should be less uncertainties in the data evaluation 
discussed in the present study. For example, in case of 
the neutron as the spectator one needs to use the miss- 
ing mass reconstruction technique. In this context let us 
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Fig. 2. Reaction cross section for (a) n^D^ n (dashed 
hne), (b) D~ p-^D~ p (sohd line) and (c) D~ p-^D'^n (dash- 
dotted line) as a function of the l)-meson momentum (lower 
axis) and the kinetic energy e in the center-of-mass system 
(cms) (upper axis). 



cesses [47, 48] . The reaction amplitude is obtained by solv- 
ing a Lippmann-Schwinger type scattering equation for 
the interaction potential. The features of the DN ampli- 
tude based on this model are much more realistic than 
the ones employed in previous studies [22-24]. Indeed, 
in the former studies the D^p cross section for instance 
was assumed to be momentum independent and equal to 
~20 mb [22,24] or 5 mb [23]. Moreover, the angular depen- 
dence of the elastic scattering was assumed to be either 
isotropic [22,24] or forward peaked [23], i.e. proportional 
to exp(6i), where t is the four momentum transfer squared, 
with a slope 6=2 GeV~^ . The reason for such assumptions 
was the lack of any microscopic calculations of the DN 
scattering amplitude in those days. 

To illustrate the differences between the previous cal- 
culations [22-24] and the results of Ref. [20] we show in 
Fig. 2 predictions for the D~ n-^D~ n, D~p D~p and 
D~p^D'^n reaction cross sections as a function of the D- 
meson momentum (lower axis) and the cm kinetic energy 
e (upper axis). It is clear that the scattering cross sec- 
tions for all channels depend significantly on the D-meson 
momentum. 

Note that the DN scattering amplitudes for the dif- 
ferent reaction channels shown in Fig. 2 are related to the 
isospin basis used in Ref. [20] by 



also recall that charged open-charm mesons are in gen- 
eral reconstructed by the leptonic, semileptonic and ha- 
dronic Kittt and Kmrn decay modes. Usually the neutral 
open-charm mesons can be well detected by their hadronic 
Z3°^A'^7r+ and D'^-^K^tt^ decays [7]. For the full recon- 
struction of the final state of the pd^D'^D~~p reaction it is 
therefore important that both semileptonic and hadronic 
modes are detected with high accuracy. This should be 
kept in mind for the design of future experiments, e.g., of 
the PANDA experiment [16]. High accuracy is crucial for 
identification of DN and DN rescattering effects, whose 
absolute values and energy dependences could be very dif- 
ferent, as indicated by several studies using different mod- 
els [22,40-44], and as will become clear from the present 
work, too. The identification of the different effects might 
be possible when a full reconstruction of the final state is 
feasible. 



3.1 The DN amplitude 

For the DN scattering amplitude we use the results of 
our recently pubhshed potential model [20]. This model 
for the DN interaction was constructed within the meson- 
exchange framework, but supplemented with a short-distance 
contribution from one-giuon-exchange. The model was de- 
veloped in close analogy to the meson- exchange KN in- 
teraction of the Jiilich group [45,46] utilizing SU(4) sym- 
metry constraints. The main ingredients of the interac- 
tion are provided by vector meson (p, lu) exchange and 
higher-order box diagrams involving D*N, DA, and D*A 
intermediate states. The short range part is supplemented 
by additional contributions from genuine quark-gluon pro- 



TM{D-n^D-n) = Tm{D"p^D"p) = /: 
TM{D-p~^D-p) 



(11) 



TMiD"n^ 



Tm{D-p- 



.D"n) = -{h 



D"n) ^ -ifo+h), (12) 
/o), (13) 



where /o and /i are the isospin /=0 and I—l amplitudes, 
respectively. 




COS tS- cos i9 

Fig. 3. Differential cross sections for the D" p^D~ p reaction 
in the cm system at different momenta. 
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Differential cross sections for D^p^D^p at different 
momenta are presented in Fig. 3. The distributions are 
almost isotropic for momenta below ~500 MeV/c, but be- 
come forward peaked at higher momenta. Note that there 
is no simple way to parametrize the angular dependence 
with functions like exp{bt) [23], unless the slope parameter 
b is taken to be momentum dependent. 

For completeness, let us mention that other models of 
the DN interaction have been published in recent years 
[41,49]. Those authors considered s- waves only. The cross 
sections predicted by these models at threshold are 8.5 mb 
[D-n _D~n), 5.54 mb {D-p D~p), and 0.03 mb 
(D-p D°n) [41,43] and 10.6 mb, 2.64 mb, and 2.64 mb 
for model B of Ref. [49] , respectively. 



3.2 The DN amplitude 

The DN interaction [21] employed in the present study is 
also constructed in close analogy to the meson-exchange 
KN model of the Jiilich group [50] as well as by exploit- 
ing the close connection between the DN and DN systems 
due to G-parity conservation. Specifically, the latter con- 
straint fixes the contributions to the direct DN interaction 
potential while the former one provides the transitions to 
and interactions in channels that can couple to the DN 
system. Accordingly, the DN interaction is likewise pro- 
vided by vector-meson (p, w) exchange and higher-order 
box diagrams involving D*N, DA, and D*A intermediate 
states. The short-ranged quark-gluon processes, however, 
are absent here because the quark-exchange mechanism 
cannot contribute to the DN interaction due to the dif- 
ferent quark structure of the D meson. As far as the cou- 
pling to other channels is concerned, we follow here the 
arguments of Ref. [50] and we take into account only the 
channels 7rylc(2285) and 7rZ'c(2455). Furthermore, we re- 
strict ourselves to vector-meson exchange and we do not 
consider any higher-order diagrams in those channels. Pole 
diagrams due to the ylc(2285) and Z'c(2455) intermediate 
states are, however, consistently included in all channels. 

In this basic model all free parameters - the coupling 
constants and the cut-off masses at the vertex form factors 
of the occurring meson-meson-meson and meson-baryon- 
baryon vertices, cf. [50] - are fixed by the assumed SU(4) 
symmetry and the connection with the KN model, re- 
spectively. When solving the coupled-channel Lippmann- 
Schwinger equation with this interaction model we observe 
that two states are generated dynamically below the DN 
threshold, one in the S'oi partial wave and the other one in 
the Sii partial wave. (We use here the standard spectro- 
scopical nomenclature Lj2j-) In view of the close analogy 
between our DN model and the corresponding KN inter- 
action [50] this is not too surprising, because also the lat- 
ter yields a quasi-bound state in the Sqi channel which is 
associated with the yl(1405) resonance. The bound states 
in both the KN and DN are generated by the strongly 
attractive interaction due to the combined effect of u, p 
and scalar-meson exchanges, which add up coherently in 
specific channels. 



It should be said that studies of the KN and DN inter- 
action within chiral unitary (and related) approaches like- 
wise generate the yl(1405) resonance dynamically but also 
states in the DN system [41,44,51]. In those approaches 
the strong attraction is also provided by vector-meson ex- 
change [41] or by the Weinberg- Tomazawa term [44,51]. 
In Refs. [42,44,49] the authors argued that the state oc- 
curing in the 5*01 channel of the charm C = 1 sector 
should be identified with the / = resonance ylc(2593). 
We adopt this viewpoint here too. Furthermore, we iden- 
tify the state we get in the Sn channel with the 1=1 
resonance Z'c(2800) [7]. 

In order to make sure that the DN model we are going 
to apply in our study of the reaction pd-^D'^D~p incor- 
porates these features also quantitatively we fine-tune the 
contributions of the scalar mesons to the DN interaction 
so that the position of those states generated by the model 
coincide with the values given in the list of the Particle 
Data Group. This can be achieved by a moderate change 
in the coupling constants of the a meson (from 1 to 2.6) 
and the ao meson (from —2.6 to —4.6), cf. Table 2 in 
Ref. [20]. 

Interestingly, our model generates a further state, name- 
ly in the Pqi partial wave, which, after the above fine- 
tuning, lies at 2803 MeV, i.e. just below the DN threshold. 
We are tempted to identfy this state with the ylc(2765) res- 
onance, whose quantum numbers are not yet established 
[7]. Though we do not reproduce the resonance energy 
quantitatively, we believe that further refinements in the 
DN model, specifically the inclusion of the AcTttt channel 
in terms of an effective aA^ channel, can provide sufficient 
additional attraction for obtaining also quantitative agree- 
ment. The mechanism could be the same as in case of the 
Roper (A^*(1440)) resonance, which is generated dynam- 
ically in the Jiilich -kN model [52,53]. Here the required 
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Fig. 4. Reaction cross sections for (a) D^n^D'^n (dashed 
line), (b) D°p-^D°p (solid line) and (c) D°p^D+n (dash- 
dotted line) as a function of the _D-meson momentum (lower 
axis) and the cms kinetic energy e (upper axis). 
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Fig. 5. Differential cross sections for the D'^p—^D^p reaction 
in the cm system at different momenta. 

strong attraction is produced via the coupling of the t:N 
p-wave (where the Roper occurs) to the s-wave in the <tN 
system, facilitated by the different parities of the tt and a 
mesons. 

Some results of our DN model are presented in Figs. 4 
and 5. The DN scattering amplitude for the different re- 
action channels shown in the Fig. 4 are related to those in 
the isospin basis by 



TM{D''n^D^n) = Tm{D+p^D+p) - A, 



(14) 



TMiD^'p^D^'p) = TMiD+n^D+n) = -(/o+/i), (15) 



TMiD"p^D+n) = -ih-fo), 



(16) 



where /o and /i are the isospin /=0 and I—l amplitudes 
respectively. 

Obviously, also the DN cross sections show a signifi- 
cant momentum dependence in all charge channels. Fur- 
thermore, the cross sections are substantially larger than 
those we obtain for DN. Specifically, for the pure 1 = 1 
channel D^n the cross section amounts to almost 600 mb 
at threshold. This is not too surprising in view of the near- 
by quasi-bound state. The latter is also reflected in the 
s-wave scattering lengths. 



7=0 
'■DN 

'■DN 



(-0.41 + i0.04) fm 
(-2.07+^0.57) fm , 



(17) 



namely by the rather large value of the real part in the 
1=1 channel. For completeness, let us mention here that 
the scattering lengths of the DN interaction of Hofmann 
and Lutz [41], reported in [43], amount to about —0.4 fm 
for both isospin channels. In agreement with that work we 



find that the imaginary part is negligibly small for / = 0. 
However, contrary to [43] in our DN model this is not the 
case for the 1=1 channel. 

Angular distributions for the reaction D^p^D^p are 
shown in Fig. 5. Obviously, there is a strong anisotropy al- 
ready at fairly low momenta. It is due to significant contri- 
butions in the Pqi partial wave in this momentum region 
induced by the near-threshold quasi-bound state produced 
by our model, as discussed above. For higher momenta the 
differential cross section becomes forward peaked, similar 
to the predictions of our model for the DN system. 

Further details of our DN model will be reported in a 
forthcoming publication [21]. 



4 DD production in pN annihilation 

The total cross section for the reaction pd — > D^D^p de- 
pends crucially on the elementary NN DD annihila- 
tion amplitude Ta- Unfortunately, so far there is no ex- 
perimental information about this reaction and even the- 
oretical studies are rather scarce [54-56]. In Ref. [54] re- 
sults were given for the reaction pp — > D^ D^ in a quark 
plus diquark model, where the elementary flavour chang- 
ing process is due to diquark-antidiquark annihilation and 
subsequent creation of a quark-antiquark pair through one 
gluon. Kaidalov and Volkovitsky calculated the reaction 
pp — > DD in the framework of a non-perturbative quark- 
gluon string model, based on secondary Regge pole ex- 
changes including absorptive corrections [55]. The result 
of both works are summarized in Fig. 6. The larger cross 
sections are predicted by the model of Ref. [54] with a 
maximal value of around 0.15 iib at piat — 12 GeV/c cor- 
responding to ^/s « 5 GeV. The model of Ref. [55] yields 
'^pp^D°D" ~ 0-05 fib at the maximum, while the one for 
pp D^D^ is a factor of 4 smaller. In both approaches, 
the unknown parameters of the models were fixed by con- 
sidering corresponding flavour changing reactions involv- 
ing strange quarks {pp KK and/or pp — > AA, etc.). 

Superficially the predictions of the two studies are qual- 
itatively similar, specifically for pp —> D^D^ which in case 
of the model of KroU et al. [54] is suppressed by roughly a 
factor four as compared to the prediction for pp — > D~ D^ 
[57] so that its maximal value would then practically co- 
incide with the one obtained in [55]. However, for the re- 
action pp D~D^ the magnitudes of the cross sections 
differ by more than a factor of 10. Furthermore, there is 
a significant difference in the energy dependence of the 
two model predictions, as seen in Fig. 6, so that the vari- 
ations of the predictions are even larger in specific energy 
regions. 

Admittedly, those results are only of limited use for 
our own investigation of the pd D~ D^p reaction. For 
example, for the amplitude corresponding to the reac- 
tion chain pd —f D^D^p — )■ D~D^p, involving rescat- 
tering in the D'^p — > D'^p and D~p —> D~p channels, one 
needs the pure 1=1 annihilation amplitude pn —> D~ D'^ 
which cannot be reconstructed from the available infor- 
mation about those models. The other contributions to 
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Fig. 6. Predictions for the pp — > DD annhilation cross section 
taken from Refs. [54] (sohd Une) and [55] (dashed and dash- 
dotted hnes). 



the reaction amplitude, involving {DN or DN) charge- 
exchange rescattering, pd — + D^D^n D^D^p and 
pd D^D^n D~D'^p, require pp D^D^ and pp 
D^D^ , respectively, but here the relative phase between 
the terms is not known. Thus, we are facing the problem 
that we either have to add all contributions incoherently 
and make additional assumptions about the isospin depen- 
dence of NN DD or we consider only the amplitude 
involving elastic DN and DN rescattering. We prefer the 
latter option. In this case we can add the Born term and 
the DN and DN rescattering contributions coherently, be- 
cause they all involve the same elementary pn D~ D^ 
annihilation amplitude, and we can include the resulting 
interference effects in the evaluation of the observables. 
However, absolute predictions are out of reach and all re- 
sults will be shown as number of events only. On the other 
hand we consider the energy dependence of the elementary 
pn — » D~ D^ annihilation amplitude in our calculation by 
adopting the results given in [55] for pp DD. But we 
should say that its influence on the observables shown in 
the present paper is practically negligible. 



5 Open charm production in pd annihilation 

Now we present results for D^D^ production in antiproton- 
deuteron annihilation utilizing the formalism and the ele- 
mentary DN and DA'' amplitudes described above, taking 
into account the Born diagram of Fig. la) and the rescat- 
tering diagram of Fig. lb). For the latter we consider both 
D~p and D^p scattering in the final state. The pn-^DD 
threshold on a free nucleon corresponds to the antiproton 
momentum of roughly 6.43 GeV/c. The absolute thresh- 
old for the DD production in antiproton-deuteron anni- 
hilation is at the antiproton momentum of 4.55 GeV/c. 
Evidently, close to the reaction threshold the production 
rate will be strongly suppressed by the phase space. We 



choose for our calculation the antiproton momentum of 
7 GeV/c, corresponding to the region where the model 
calculation of [55] predicts the largest cross sections for 
the elementary NN DD reaction. 



5.1 Spectator momentum distribution 

In Fig. 7 we present our predictions for the spectator 
proton momentum distribution. Here the solid histogram 
indicates the full result that includes the Born and the 
rescattering diagrams while the dashed line is the result 
based on the nucleon-exchange Born diagram alone, both 
obtained with the deuteron wave function of the CD Bonn 
NN potential [37] . Since the absolute normalization of the 
reaction cross section is quite uncertain we show the re- 
sults as number of generated events. 

The dotted and dash-dotted curves are results for the 
Born term alone employing the deuteron wave functions 
of the Paris [38] and full Bonn [39] potentials, respectively. 
Obviously, there is some model dependence which becomes 
more pronounced for spectator proton momenta above 300 
MeV/c. But the rescattering mechanism is definitely by 
far the most dominant effect for momenta from around 400 
MeV/c upwards. Since rescattering occurs in both D^p 
and D^p systems one needs to apply specific methods to 
separate their contributions in a reliable way, as will be 
discussed below. 

In comparison to the multipion production case, cf. 
Figs. 11 and 12 in the Appendix, the enhancement due to 
the rescattering processes sets in at noticably higher spec- 
tator momenta and is also less pronounced. It was argued 
[25,26] that the strong enhancement seen in the proton 




0.8 1 
p, (GeV/c) 

Fig. 7. Proton momentum spectrum for the pd~*D" D" p re- 
action. The dashed (dotted, dash-dotted) line shows the result 
of a calculation for the nucleon exchange Born diagram only 
(Eq. (2)) based on the s and d-wave parts of the deuteron wave 
function of the CD Bonn [37] (Paris [38], full Bonn [39]) NN 
potential. The solid histogram is the full calculation (for CD 
Bonn) that includes D'^p and D~p rescattering. 
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Fig. 8. Dahtz plot for the pd^D^' p reaction. Here, M{Dp) and M{DD) are the D^p and D'^ invariant masses, re- 
spectively. The arrows points in the direction of decreasing intensity of the distribution. Each plot a)-d) is explained in the 
text. 



momentum spectrum for the multipion reactions is to a 
good part due to the excitation of the Z\(1232) resonance 
in the ttN rescattering processes. Though our DN scat- 
tering amplitude is dominated likewise by poles, in several 
partial waves, cf. the discussion above, their influence on 
the momentum spectrum seems to be smaller, presumably 
because they all lie below the DN elastic threshold. The 
enhancement we get for the pd^D^D~p reaction seems to 
be somewhat smaller than what was reported in an earlier 
model calculation by Gassing et al. [23]. But one has to 
keep in mind that in the latter work the contribution of 
the cZ-wave component to the Born (spectator) term was 
neglected and, moreover, the results for the Born term 
and the rescattering term are shown separately, while we 
added them coherently in our calculation. 



5.2 Dalitz plot 

A well-known method for the reconstruction of the re- 
action dynamics is the Dalitz plot analysis of the final 
state [33,58-61]. For instance, for multi-pion production 
from pd annihilation the presence of rescattering effects 
was demonstrated via the projection of the Dalitz plot on 
the invariant mass spectrum of the final 7rp system [32, 
62] . An analysis in form of a partial wave decomposition 
of the Dalitz plot was proposed [58-60] for the pp-^pK^ A 
reaction, which finally allowed to study and separate [61] 
non-resonant and resonant contributions in the K~^A sub- 
system. A similar technique could be applied in the anal- 
ysis of the pd^D'^D~p reaction. This method allows to 
study all subsystems, D°p, D~p and D^D~ , but obviously 
requires high mass resolution and significant statistics. 

Fig. 8 presents the Dalitz plot evaluated for the reac- 
tion pd-^D^D~p at the antiproton momentum of 7 GeV/c. 
The horizontal axes in the plot indicate the invariant mass 



of the D'^D^ system, while the vertical axes indicate the 
mass of the D~p system. Here panel a) shows the results 
obtained with the nucleon-exchange Born diagram only, 
b) those obtained with D~p rescattering alone, while c) 
illustrates the results with D^p FSI alone. Finally, Fig. 8d) 
contains the full results, i.e. when all three contributions 
are included coherently. The arrows in the figures indi- 
cate the direction of decreasing intensity of the distribu- 
tion. Note that we implemented a cut on the spectator 
proton momentum in the calculations in order to reduce 
the contribution from the Born diagram. Specifically, we 
considered only events involving spectator protons with 
momenta above 300 MeV/c. 

Evidently, the differences between the distributions re- 
sulting from the different diagrams are quite significant. 
The result based on the Born term alone indicates strong 
correlations between the D^- and D~-meson invariant mass 
This is due to the fact that they are produced from the 
same vertex, namely via. pn—^ D~ . The invariant energy 
of the D^D~ system is essentially given by the energy of 
the incoming antiproton, while the dispersion of the dis- 
tribution is related to the square of the deuteron wave 
function. Assuming the target neutron to be at rest, the 
invariant mass of the D^D~ system produced in the re- 
action pn—>D'^D~ at antiproton momentum of 7 GeV/c 
is equal to 3.86 GeV. The strong D'^D^ correlation pro- 
duces also a kinematical reflection, detectable in the D~p 
system, namely in form of an enhancement in the high 
mass D^p spectrum. However, this enhancement, being 
purely kinematical, does not have any relevance for the 
interpretation of the rescattering mechanism. 

Results obtained for the rescattering diagrams alone 
are shown in Figs. 8b) and c). It is clear that there are no 
strong correlations in the D^D~ system anymore. Now 
the spectrum is primarily distorted by the corresponding 
rescattering terms. The D~p and D'^p projections of the 
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Fig. 9. The missing mass distribution of Px given by Eq. (18) 
for the pd~*D° D^p reaction. The dashed histogram shows re- 
suhs obtained with the nucleon-exchange Born diagram and 
the D^p FSI for spectator proton momenta above 300 MeV/c. 
The sohd histogram is the result for the D~p rescattering di- 
agram based on Eq. (3). 

Dalitz plot show the distribution produced by the relevant 
scattering amplitudes. 

The final distribution, shown in Fig. 8d), corresponds 
to a calculation that includes the Born diagram plus both 
rescattering diagrams. It is fairly non-uniform and clearly 
indicates substructures resulting from the individual reac- 
tion mechanisms. 



5.3 Missing mass of the exchanged meson 

A method [63] that could be useful for the separation of 
the DN and DN rescattering contributions is based on 
the assumption that the dominant part of the rescattering 
amplitude comes from contributions where the particles in 
the intermediate state are on shell. We should emphasize 
that this method has no physical meaning when it comes 
to the off-shell part of the rescattering diagrams because 
then the exchanged meson is virtual. 

Under the assumption that the charmed meson and 
nucleon are on-shell before undergoing rescattering, one 
can reconstruct the four-momentum of the meson in the 
loop following the missing mass technique via 

Pi = (p, + Pz)- - ^^7v)^ (18) 

where Pg and Pq- are the four- momenta of the specta- 
tor proton and the final meson, and P/v is the four- 
momentum of the nucleon in the loop (i.e. the one involved 
in the rescattering process) , given by the loop momentum 
q and energy Eg as in Eq. (3). Let us assume for the mo- 
ment that the scattering process takes place on a free nu- 
cleon at rest. Then we would get P^=mjj, where mjj is 



the mass of the (incoming) D meson. However, since the 
reaction does not take place on a free nucleon but on a 
nucleon from the initial deuteron, the interacting nucleon 
is not at rest and, therefore, one expects a distribution 
of the missing mass in Eq. (18) around the central value 
oi tud reflecting the Fermi motion of the nucleon in the 
deuteron. 

It is clear that the Born term as well as D^p scatter- 
ing would not lead to such a distribution because in this 
case there is no correlation between the four-momenta in 
Eq. (18). 

In Fig. 9 we present results for the missing mass of the 
exchanged meson as given by Eq. (18), obtained for the 
antiproton momentum of 7 GeV/c. In the corresponding 
calculations the Born diagram and the rescattering ampli- 
tudes according to Eq. (3) are taken into account. Then 
we evaluate Eq. (18) using the four-momenta of the final 
spectator proton and of the final Z?^ -meson, and assume 
that PN — {mp, 0). The dashed histogram in Fig. 9 includes 
contributions from the nucleon-exchange Born diagram for 
spectator momenta above 300 MeV/c as well as from the 
D°p FSI. The solid histogram is the result obtained for 
the D~p rescattering term. 

The results shown in the Fig. 9 look very promising 
with regard to the possibility for a separation of the re- 
action mechanisms. However, one should keep in mind 
that there are uncertainties due to the on-shell assump- 
tion made in the model calculation as well as in the eval- 
uation of the missing mass, which cannot be quantified 
easily. Nonetheless, we want to mention that this method 
was actually used in Ref. [63] in the data evaluation and 
reconstruction of the hyperon production mechanisms in 
antiproton annihilation on xenon nuclei. Recently, this 
missing-mass method was also utilized for the analysis of 
K^K~ pair production from carbon [64]. It is important 
to stress that the method cannot and should not be ap- 
plied for too high momenta of the spectator proton, where 
the reaction is dominated by off-shell contributions and, 
therefore, the basic assumptions of the method are evi- 
dently no longer valid. 

We have also performed calculations of the missing 
mass distribution for the corresponding case of a final D'^ 
meson where one can isolate the D'^p rescattering contri- 
butions. The results are qualitatively very similar to the 
D~ case and therefore we do not show them here. 



5.4 Correlation between the scattering planes 

Finally, we discuss the correlation between the two scat- 
tering planes [65]. One plane is given by the momenta 
of the antiproton and of the spectator proton. The other 
one is fixed by the momenta of the antiproton and of 
the produced charmed meson, the D~ meson, say. Then, 
due to the conservation of the transverse momenta in the 
D~p^D~p scattering process the azimuthal angle be- 
tween these planes is peaked around 0~18O°. This cor- 
relation simply follows from the reaction kinematics. If 
the spectator proton is at rest before the scattering and 
the Z?~-meson has no transverse momentum, then, after 
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Fig. 10. The distribution of the azimuthal angle between the 
scattering plane given by the momenta of the antiproton and 
spectator proton and the plane fixed by the momenta of the an- 
tiproton and the D~ -meson. The results are for spectator pro- 
tons with momenta above 300 MeV/c. The upper panel shows 
the result for the Born diagram alone, while in the lower panel 
calculations including the D~p and D"p rescattering diagrams 
are presented. Here the solid histogram shows the distribution 
obtained for rescattering of the D~ -meson, while the dashed 
histogram is the corresponding result for rescattering of the 
_D'^-meson. 



D^p—^D^p scattering the transverse components of the 
momenta of the final Z)~ -meson and proton must be ex- 
actly the same but aligned in opposite direction. However, 
both the Fermi motion in the deuteron and the pn^D^ 
annihilation allow for some variations in the transverse 
momenta of the spectator proton and of the D~-meson. 
That is why in an actual experiment one would expect 
a distribution of the azimuthal angle around the value 
0=180°. 

Results of our model calculation for the distribution 
of the azimuthal angle are presented in Fig. 10. They are 
obtained again by imposing a cut on the spectator pro- 
ton momentum so that only momenta above 300 MeV/c 
contribute. The scattering plane is fixed by the momenta 
of the antiproton and the Z?~-meson. The upper panel of 
Fig. 10 shows predictions for the nucleon-exchange Born 
diagram while the lower panel contains results including 
the rescattering diagrams. Here the solid histogram corre- 
sponds to the contributions of the D^p rescattering dia- 
gram and the dashed histogram to those from D^p rescat- 
tering. From these results it seems feasible that the two 
rescattering contributions can be well isolated. 

Also here we have considered the corresponding case 
for the D'^ meson, where again the results turned out to 
be qualitatively similar. 



6 Summary 

In this paper we examined the possibility to extract in- 
formation about the DN and DN interactions from the 
pd-^D'^D~p reaction. We utilized the notion that those 
open-charm mesons are first produced by annihilating an- 
tiprotons on the deuteron and subsequently rescatter on 
the remaining (spectator) nucleon. The latter process is 
then exploited for investigating the DN and DN interac- 
tions. To explore the potential of a corresponding exper- 
iment we performed concrete model calculations taking 
into account the nucleon-exchange Born diagram as well 
as rescattering diagrams. 

As a test of the approach we first applied it to mul- 
tipion production in pd annihilation and we compared 
our results with data available for the pd— !-7r+27r~p and 
pd^2TT'^3TT~p reactions. These data [66-68] show strong 
evidence for contributions from ttN—htN rescattering, which 
can be seen in the spectra of the spectator proton and in- 
variant mass distribution of the final pion and proton. We 
obtained very reasonable agreement with the data avail- 
able for the spectator proton distributions. 

In our investigation of the pd-^D^D~p reaction we 
utihzed reahstic DN^DN [20] and DN^DN [21] scat- 
tering amplitudes. This is a substantial improvement over 
previous studies [22-24] which employed simplistic DN 
and DN scattering amplitudes based on somewhat ques- 
tionable assumptions. 

We found that below spectator momenta of around 
300 MeV/c the reaction is dominated by the nucleon- 
exchange Born diagram. For higher spectator momenta 
there is a sizable contribution from the rescattering di- 
agrams. In particular, their contribution is significantly 
larger than the uncertainties due variations in the high- 
momentum component of the deuteron wave function. Thus, 
selecting events with spectator momenta above 300 or 
400 MeV/c, say, should allow to obtain a data sample that 
can be used for extracting information about the DN and 
DN interactions. 

Subsequently we explored different methods for isolat- 
ing the contributions from the DN and DN rescatter- 
ing terms. We showed that the missing mass technique 
and the correlation between the planes given by the scat- 
tered meson and nucleon allow a reasonable reconstruction 
of the reaction dynamics and to separate the contribu- 
tions of DN rescattering from those of DN rescattering. 
Since these methods are based on the reaction kinematics 
we consider them as promising tools to extract informa- 
tion on the DN and DN interactions from the reaction 
pd^D"D-p. 
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A Multi-pion production in pd annihilation 

To illustrate the applicability of the discussed formalism 
we consider experimental data available for the reactions 
pd— >27r"''37r~p and pd— i-7r+27r^p obtained at the Low En- 
ergy Antiproton Ring (LEAR) at CERN using annihila- 
tion at rest in a hydrogen gas [66, 67] and at the Brookhaven 
National Laboratory (BNL) using the deuterium bubble 
chamber [68]. 

We calculate the proton spectator distribution by sum- 
ming the Born and rescattering amplitudes and integrat- 
ing over the 6-body and 4-body final states. The calcu- 
lations are done with the deuteron wave function of the 
CD-Bonn potential [37]. Available data [69,70] on pion 
multiplicities for pp annihilation at rest and in flight show 
practically no dependence on the antiproton momentum 
within the range up to ~100 MeV/c. Therefore, we as- 
sume the annihilation amplitude Ta to be a constant. The 
spectator proton momentum distribution was measured in 
different experiments [66-68] and the data were published 
with arbitrary normalization depending on the total num- 
ber of detected events. Therefore, we normalize our calcu- 
lation to the data but we also normalize the different data 
sets to each other, as explained below. 

For evaluating the contribution from pion-nucleon re- 
scattering we use the current solution of the GWU/CNS 
partial wave analysis [71,72]. We account for isospin and 
topological factors following the prescription given in 
Refs. [30,31,35]. The phase-space integration is done by 
the Monte-Carlo method based on event-by-event simu- 
lations. This allows us to apply kinematical cuts on the 
final pion momenta similar to that discussed in Refs. [66, 
67] in order to investigate discrepancies between the data 
at spectator proton momenta below 250 MeV/c. We will 
come back to this issue later. 

Experimental results for the proton momentum distri- 
bution for the reaction pd—!-27r"*'37r~p are shown in Fig. 11. 
The squares [66] and circles [67] are data from the LEAR 
facility while the triangles are from an experiment [68] 
at the Brookhaven National Laboratory (BNL) using the 
deuterium bubble chamber. 

The basic difficulty in the counter measurements [66, 
67] is the reconstruction of the low-momentum part of the 
spectator spectra. As mentioned in Ref. [66], for the di- 
rect proton detection at least transverse proton momenta 
of 130 MeV/c are required. The very low momentum pro- 
tons can be only reconstructed through exclusive measure- 
ments of the final pions and applying the missing momen- 
tum method. But such a reconstruction introduces addi- 
tional uncertainties in the low momentum spectator spec- 
tra. To avoid any ambiguity we normalized all data sets 
at momenta around 400 MeV/c, i.e. at values where the 
proton spectator momentum was measured directly. This 



normalization emphasizes that the shape of the measured 
spectrum from the different experiments is almost iden- 
tical at higher spectator proton momenta. On the other 
hand, we see a substantial disagreement between the avail- 
able data [66-68] at momenta below ~ 250 MeV/c. 

Experimental results for the reaction pd— >7r+27r^p are 
presented in Fig. 12. Again, the squares are LEAR data [66] 
while the triangles are from the BNL [68] . Both data sets 
are normalized at proton momenta around 400 MeV/c 
and in such a way that the scale is roughly the same as in 
Fig. 11. This allows us to compare the shapes of the spec- 
tator proton distributions for the two reactions. There is 
clearly a difference between those shapes for pd— >27r+37r~p 
and pd— !-7r+27r~p. Indeed one expects that the shape of the 
spectator momentum distribution depends on the momen- 
tum carried by the scattering meson and, consequently, 
that reactions with different final pion multiplicity would 
exhibit different shapes of the proton spectra. The more 
energetic pions from the pd— >7r+27r~p reaction produce 
more energetic spectator protons in the rescattering. Qual- 
itatively this follows from the rescattering amplitude of 
Eq. (3).^ 

Again, for pd— >7r+27r p the experimental results are 
consistent at proton momenta above 250 MeV/c and dis- 
agree substantially at low momenta. As already indicated 
above, we applied kinematical cuts on the final pion mo- 
menta, similar to those discussed in Refs. [66,67], in the 
course of our investigation for the reaction pd—>7r+27r~p as 
well as forpd^27r+37r~p. But it turned out that those cuts 
do not resolve the disagreement between the data [66-68] 
at spectator proton momenta below 250 MeV/c. 




Ps (GeV/c) 



Fig. 11. Proton momentum spectrum for the pd— >27r+37r~p re- 
action. The data are from Refs. [66] (squares), [67] (circles) and 
[68] (triangles). The dashed line shows the result from Eq. (2) 
taking into account s- and d-wave parts of the deuteron wave 
function while the dotted line is based on the s-wave compo- 
nent alone. The solid histogram is the full calculation including 
the nucleon exchange Born diagram and rescattering diagrams 
using the GWU/CNS TriV partial- wave amplitudes [71,72] for 
Tm. 
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Fig. 12. Proton momentum spectrum for the pd-^-rt^ 2tx^ p 
reaction. The data are from Ref. [66] (squares) and [68] (trian- 
gles). Same description of curves as in Fig. 11. 

Let us now come to the results of our model calculation 
and first discuss the normalization, which is a somewhat 
delicate issue. From a theoretical point of view it should 
be done preferrably around the peak in the distribution 
at low momenta where the spectrum is dominated by the 
Born diagram and the s-wave component of the deuteron 
wave function. But this is exactly the region where the 
experimental uncertainty is very large. Thus, we decided 
to normalize our results also in the plateau region, i.e. 
around 400 MeV/c. Note that the normalization is done 
for the full calculation. The relative size of the momen- 
tum distribution at the peak as compared to the plateau 
is fixed by the ingredients of the model alone, i.e. the con- 
tributions of the Born diagram and from the rescattering 
diagram. There is no additional normalization constant 
involved here. 

The dashed lines in Figs. 11 and 12 correspond to 
the contribution from the nucleon-exchange Born diagram 
given by Eq. (2). Interestingly, our predictions for low mo- 
menta agree well with the data of [67] (for pd-^27r^3Tr^p) 
and [66] (for pd—nT^2'K~p). For proton momenta above 
200 MeV/c all data show a substantial enhancement with 
respect to the predictions based on the Born term alone. 
For illustration purposes we show here also results using 
only the s-wave part of the deuteron wave function (dot- 
ted line). There have been some speculations that the en- 
hancement at higher momenta could indicate an excitation 
of the short range component of the deuteron wave func- 
tion [73-77]. But, in any case, a simple renormalization 
of the d-wave contribution would not reproduce the ob- 
served shape of the spectator proton distribution around 
the plateau, i.e. for ps~400 MeV/c. 

The solid histograms in Figs. 11 and 12 are the re- 
sults of our full calculation including the nucleon-exchange 
Born diagram and the rescattering diagram. But we should 
say that the off-shell corrections of Eq. (6), not considered 
in the present investigation, are known to lead to varia- 



tions of the order of 30% or more in the absolute value 
of the rescattering contribution, though they do not effect 
the shape of the spectator proton distribution [31]. 

Our results are in reasonable agreement with that of 
Ref. [31], but are in contradiction to the conclusions of 
Ref. [78], where the spectator proton momentum distri- 
bution is well reproduced by taking into account the Born 
diagram of Fig. la) and pion absorption on the spectator 
nucleon. However, there is strong experimental evidence in 
favor of the rescattering mechanism: The invariant mass 
spectrum of the 27r+27r~ system from the pd— !-27r+37r~p 
reaction was measured [67] for different cuts on the spec- 
tator proton momenta. If the final pions do not undergo 
rescattering, such cuts should not change the invariant 
mass distribution. But in the experiment it turned out 
that the 27r+27r^ invariant mass distribution depends sub- 
stantially on the proton momentum cut when taken below 
or above Ps = 200 MeV/c. 

References 

1. T. Matsui and H. Satz, Phys. Lett. B 178, 416 (1986). 

2. H. Satz, J. Phys. G 32, R25 (2006) [arXiv:hep- 
ph/0512217]. 

3. S. Batsouli, S. Kelly, M. Gyulassy and J.L. Nagle, Phys. 
Lett. B 557, 26 (2003) [arXiv:nucl-th/0212068]. 

4. M. Djordjevic and M. Gyulassy, Phys. Lett. B 560, 37 

(2003) [arXiv;nucl-th/0302069]. 

5. N. Armesto, C.A. Salgado and U.A. Wiedemann, Phys. 
Rev. D 69, 114003 (2004) [arXiv:hep-ph/0312106]. 

6. M. Djordjevic and M. Gyulassy, Nucl. Phys. A 733, 265 

(2004) [arXiv:nucl-th/0310076]. 

7. W.-M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 
(2006) 

8. G.D. Moore and D. Teaney, Phys. Rev. C 71, 064904 

(2005) [arXiv:hep-ph/0412346]. 

9. H. van Hees and R. Rapp, Phys. Rev. C 71, 034907 (2005) 
[arXiv:nucl-th/0412015] . 

10. H. van Hees, V. Greco and R. Rapp, Phys. Rev. C 73, 
034913 (2006) [arXiv:nucl-th/0508055 

11. S. Wicks, M. Djordjevic, C. Horowitz and M. Gyulassy, 
Nucl. Phys. A 784, 426 (2007) [arXiv:nucl-th/0512076]. 

12. S.S.Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 
94, 082301 (2005)[arXiv:nucl-ex/0409028]. 

13. J. Bielcik et al. [STAR Collaboration], Nucl. Phys. A 774, 
697 (2006) [arXiv:nucl-ex/0511005]. 

14. Y. Zhang et al. [STAR Collaboration], J. Phys. G 32, S529 

(2006) [arXiv:nucl-ex/0607011]. 

15. C. Zhong et al. [STAR Collaboration], J. Phys. G 34, S741 

(2007) [arXiv;nucl-ex/0702014]. 

16. M. KotuUa et al.. Technical Progress Report for PANDA, 
Strong Interaction Studies with Antiprotons, Febru- 
ary 2005, http://www-panda.gsi.de/db/papersDB/PC19- 
050217_panda_tpr.pdf. 

17. O.N. Hartmann et al. [PANDA Collaboration], Int. J. 
Mod. Phys. A 22, 578 (2007). 

18. K.T. Brinkmann, P. Gianotti, I. Lehmann [PANDA 
Collaboration], Nucl. Phys. News 16, 15 (2006) 
[arXiv:physics/0701090]. 

19. W. Kuhn [PANDA Collaboration], Acta Phys. Polon. B 
37, 129 (2006). 



J. Haidenbauer et al.: Charmed meson rescattering in the reaction pd — > DDN 



13 



20. J. Haidenbauer, G. Krcin, Ulf-G. Meifiner and A. Sibirtsev, 
Eur. Phys. J. A 33, 107 (2007) [arXiv:nucl-th/0704.3668]. 

21. J. Haidenbauer et al., in preparation. 

22. A. Sibirtsev, K. Tsushima aud A. W. Thomas, Eur. Phys. 
J. A 6, 351 (1999) [arXiv:uucl-th/9904016]. 

23. W. Gassing, Y. S. Golubcva and L. A. Kondratyuk, Eur. 
Phys. J. A 7, 279 (2000) [arXiv:nucl-th/9911026l. 

24. A. Sibirtsev, Nucl. Phys. A 680, 274c (2001). 

25. C.G. Fasano and M.P. Locher, Z. Phys. A 336, 469 (1990). 

26. V.M. Kolybasov, I.S. Shapiro and Yu.N. Sokolskikh, Phys. 
Lett. B 222, 135 (1989). 

27. J.M. Laget, Phys. Kept. 69, 1 (1981). 

28. G.G. Fasano and M.P. Locher, Z. Phys. A 338, 197 (1991). 

29. J.M. Laget, Nucl. Phys. A 296, 388 (1978). 

30. M.P. Locher and B.S. Zou, Z. Phys. A 340, 187 (1991). 

31. G.G. Fasano, M.P. Locher and S. Nozawa, Z. Phys. A 338, 
95 (1991). 

32. D.V. Voronov and V.M. Kolybasov, JETP Lett. 57, 162 
(1993). 

33. E. Byckling and K. Kajantie, Particle Kinematics (Wiley 

and Sons, New York, 1973). 

34. J.M. Laget, Phys. Rev. G 73, 044003 (2006). 

35. M.P. Locher and B.S. Zou, J. Phys. G 19, 463 (1993). 

36. J.M. Laget, Phys. Rev. G 75, 014002 (2007) [arXivinucl- 
th/0603009]. 

37. R. Machlcidt, Phys. Rev. G 63, 024001 (2001) [arXivmucl- 

th/0006014]. 

38. M. Lacombe, B. Loisoau, J. M. Richard, R. Vinh Mau, 
J. Gote, P. Pires and R. De Tourreil, Phys. Rev. G 21, 861 
(1980). 

39. R. Machleidt, K. Holinde, and Gh. Elster, Phys. Rept. 149, 

1 (1987). 

40. Z. w. Lin, G. M. Ko and B. Zhang, Phys. Rev. G 61, 024904 
(2000) [arXiv:nucl-th/9905003]. 

41. J. Hofmann and M. F. M. Lutz, Nucl. Phys. A 763, 90 

(2005) [arXiv:hep-ph/0507071]. 

42. L. Tolos, J. Schaffner-Bielich and A. Mishra, Phys. Rev. G 
70, 025203 (2004) [arXiv:nucl-th/0404064]. 

43. M. F. M. Lutz and G. L. Korpa, Phys. Lett. B 633, 43 

(2006) [arXiv:nucl-th/0510006]. 

44. T. Mizutani and A. Ramos, Phys. Rev. G 74, 065201 
(2006) [arXiv:hep-ph/0607257]. 

45. R. Biittgon, K. Hohndc, A. Miillcr-Grochng, J. Speth, and 
P. Wyborny, Nucl. Phys. A 506, 586 (1990). 

46. M. HofTmann, J.W. Durso, K. Hohnde, B.G. Pearce, and 
J. Speth, Nucl. Phys. A 593, 341 (1995). 

47. D. Hadjimichef, J. Haidenbauer and G. Krein, Phys. Rev. 
G 66, 055214 (2002) [arXiv:micl-th/0209026]. 

48. J. Haidenbauer and G. Krcin, Phys. Rev. C 68, 052201 
(2003) [arXiv:hcp-ph/0309243]. 

49. L. Tolos, A. Ramos and T. Mizutani, Phys. Rev. G 77, 
015207 (2008) [arXiv:0710.2684 [nucl-th]]. 

50. A. Miiller-Groehling, K. Holinde, and J. Speth, Nucl. Phys. 
A 513, 557 (1990). 

51. M. F. M. Lutz and E. E. Kolomeitsev, Nucl. Phys. A 730, 
110 (2004) [arXiv:hep-ph/0307233]. 

52. O. Krehl, G. Hanhart, S. Krewald and J. Speth, Phys. Rev. 
G 62, 025207 (2000) [arXiv:nucl-th/9911080]. 

53. A. M. Gasparyan, J. Haidenbauer, G. Hanhart and 
J. Speth, Phys. Rev. G 68, 045207 (2003) [arXivmucl- 
th/0307072]. 

54. P. KroU, B. Quadder and W. Schweiger, Nucl. Phys. B 
316, 373 (1989). 



55. A. B. Kaidalov and P. E. Volkovitsky, Z. Phys. G 63, 517 

(1994) . 

56. B. Kerbikov and D. Kharzeev, Phys. Rev. D 51, 6103 

(1995) [arXiv:hep-ph/9408378]. 

57. W. Schweiger, private communication. 

58. G. Hanhart, Phys. Rept. 397, 155 (2004) [arXivrhep- 
ph/0311341]. 

59. A. Sibirtsev, J. Haidenbauer, H.-W. Hammer and S. 
Krewald, Eur. Phys. J. A 27, 269 (2006) [arXivinucl- 

th/0512059]. 

60. A. Sibirtsev, J. Haidenbauer, U.-G. Meifiner Phys. Rev. 
Lett. 98, 039101 (2007) [arXiv:hep-ph/0607212]. 

61. S. Abdel-Samad et al. [GOSY-TOF Collaboration] Phys. 
Lett. B 632, 27 (2006). 

62. A. Sibirtsev, K. Tsushima and A. Faessler, Z. Phys. A 354, 
215 (1996) [arXiv:nucl-th/9511008]. 

63. V.V. Barmin et al. [DIANA GoUaboration] , Nucl. Phys. A 
683, 305 (2001). 

64. M. Hartmann, Yu. Kiselev et al. [ANKE Collaboration], 
COSY Proposal 147. 

65. Y.S. Golubeva, W. Gassing, L.A. Kondratyuk, A. Sibirtsev 
and M. Biischer, Eur. Phys. J. A 7, (2000) [arXiv:nucl- 
th/9910028]. 

66. J. Riedlberger et al., Phys. Rev. G 40, 2717 (1989). 

67. S. Ahmad et al, IV*'' LEAR Workshop: Physics at LEAR 

with Low-Energy Antiprotons, Villars-sur-OUon, Switzer- 
land (1987), edited by G. Amsler et al., Nucl. Sci. Res. 
Gonf. Ser. 14, 447 (1988). 

68. P.D. Zemany, Z.M. Ma and J.M. Mountz, Phys. Rev. Lett. 
38, 1443 (1977). 

69. J. Vandcrmculcn, Z. Phys. C 37, 563 (1988). 

70. G. Amsler, Rev. Mod. Phys. 70, 1293 (1998) [arXiv:hcp- 
ex/9708025]. 

71. R.A. Arndt, I.I. Strakovsky and R.L. Workman Phys. Rev. 
G 52, 2120 (1995) [arXiv:nucl-th/9505040]. 

72. R.A. Arndt, W.J. Briscoe, I.I. Strakovsky and R.L. 
Workman, Phys. Rev. G 74, 045205 (2006) [arXivinucl- 

th/0605082]. 

73. P. Benz and P. Soding, Phys. Lett. B 52, 367 (1974). 

74. P.J. Mulders and A.W. Thomas, Phys. Rev. Lett. 52, 1199 
(1984). 

75. Y.E. Kim and M. Orlowski, Phys. Rev. C 29, 2299 (1984). 

76. L.Ya. Glozman, Prog. Part. Nucl. Phys. 34 123 (1995). 

77. L.-C.Lu and T.-S. Cheng, Phys. Lett. B 386, 69 (1996). 

78. A.E. Kudryavtsev and V.E. Tarasov, Sov. J. Nucl. Phys. 
54, 36 (1991). 



